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1 . 


INTRODUCTION 


In this paper we study the global qualitative 
behavior of axisymmetric buckled states of homogeneous 
isotropic nonlinearly elastic shells that can suffer 
flexure, compression, and shear. Our model is geometrically 
exact in the sense that a geometric quantity, such as sm0, 
is not replaced by an approximation, such as 0 or 0 - 0V6. (The 
usual justification for such a replacement is that 0 is known 
to be small for the physical situation under study. But it 
is quite possible that a mathematical model with exact 
geometry permits only small solutions, while those with approx- 
imate geometry have large solutions.) We allow the material 
properties to be described by a very general class of non- 
linear constitutive relations. Consequently our governing 
equations form a quasilmear sixth-order system of ordinary 
differential equations. 

Our main result is that nontrivial branches of solutions 
pairs can be globally characterized by a novel nodal pattern 
relating shear and bending effects. The detection of this 
pattern is greatly aided by the generality of our theory, 
which does not obscure the simple mathematical structure of 
the governing equations with approximate geometrical relations. 

indeed, the beautiful numerical results of Bauer, Reiss, 

& Keller (L970) for their technical theory of shells indicate 
that their solution branches do not enjoy the nodal properties 
we discover for our exact theory. We surmise that no such 
technical theory would have these properties. Consequently, 
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the nature of the bifurcation diaqram away from the trivial 
solution for these technical theories could differ markedly 
from ours. 

There is a large literature on the buckling of spherical 
shells. The articles by Bauer, Reiss, & Keller (1970) , 

Antman (1971), Lange & Kriegsmann (1981), and shih (1985) 
give extensive lists of references. Our results differ 
markedly from those developed in these references (except 
for the work of Antman) by virtue of the generality of our 
theory and the global character of our theorems. 

The heart of our paper is m Sections 5 and 6. In 
Section 5 we show that our equations meet the hypotheses of 
Rabinowitz's Global Bifurcation Theorem. In Section 6 we 
prove a basic uniqueness theorem that enables us to show 
that solution branches preserve the nodal pattern they inherit 
from the trivial branch. The use of the isotropy of the 
shell material is crucial m enabling us to overcome the 
difficulties posed by the polar singularities. 

Notation . Partial derivatives are denoted by subscripts. 
Ordinary derivatives with respect to the variable s are 
denoted by primes. If f and g are functions of x and y, then 

d - denotes the m atrix of partial derivatives of f and g 
3 (x,y) 


with respect to x and y. 
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§2. FORMULATION OF THE GOVERNING EQUATIONS . 

Let { i,j,k } be a fixed right-handed orthonormal 
basis for the Euclidean 3-space E . For each <j> in R 
we define 

(2.1) ei (4>) = cos <j>i + sin 4>j, e 2 (<j>) = -sin <j>i + cos 4>j, 
ea = k 

To each (s,<J>) e [0 ,tt] x [0,2tt] there corresponds exactly 
one point on the sphere of radius 1 centered at the origin 
with position vector of the form 

(2.2) r*(s,<f>) = sin s e^cfc) - cos s e 3 

Note that s measures the arc length to £*(s,<j>) from the 
south pole of the sphere. This convention will simplify 
some of the formulas. We interpret the sphere defined by 

(2.2) to be the natural reference state of the midsurface 
of a spherical shell. The coordinates (s,4>) identify 
material points on this surface. 

An axisymmetric -configuration of a spherical shell 
that can suffer flexure, extension, and shear is determined 
by a pair of vector functions 

(2.3) (s,<j>) h- r (s, <(>) , b (s, <f>) 

with r(s,*) and b(s,*) having period 2tt and with 
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(2.4) r (s, 4>) • e 2 (4>) = 0 , 
b(s,<j>) • e 2 (<f>) =0. 

| b (s, <J>) |=1 

The reference configuration of the shell is given by 
r = r* , b = -r* . The vector b(s,<f>) is interpreted as 
characterizing the deformeff configuration of the material 
points that in the natural state of the shell lie along the 
ray determined by r*(s,<J>) 

We set 


(2.5) 

r(s,<j>) = 

r(s)ei (<j>) + z (s)k , 

(2.6) 

a(s,4>) = 

cos \p (s) ei (4>) + sin i|Ms)k , 


b (s, 4>) = 

- sin ^(s)ei(<j>) + cos iMs)k , 

(2.7) 

r s (s,<j>) = 

v (s) a (s , <j>) + n (s)b (s,4>) , 

(2.8) 

T (s) = 

r (s) 

... • _ g 



sin s 


y (s) = 

ip' (s) - 1 , 


a (s) = 

sin ( s ) ^ 



sin s 

The strain 

variables 

for our problem are 

(2.9) 

w = ( V, 

x, n, y, a) 

Let 

ni ( s , 4> ) 

and mi(s,$) denote the resultant 


contact force and contact couple per unit reference length 
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of the circle <j> r*(s,cj>) of radius sin s that are 

exerted across this circular section at r*(s,<J>) . Let 

n 2 (s,<j)) and m 2 (s,<j>) denote the resultant contact force 
and contact couple per unit reference length of the circle 
s w r*(s,<j>) of radius 1 that are exerted across this 
section at r*(s,<j>) . Since we seek axisymmetric solutions, 
we require these resultants to have the form 


(2.10) ni(s,<f>) = N (s) a (s, <(>) + H(s)b(s,<J>) , 

A 

(2.11) n 2 ( s , ) = T(s)e 2 (<M / 

A 

(2.12) mi(s,4>) = - M(s)e 2 (4>) , 

A y\ 

(2.13) m 2 (s,<})) = Z(s)cos ( s ) e x ( ) + L(s)k 


Then the equilibrium equations for the shell under an external 
hydrostatic pressure of intensity A per unit actual area 
are 


(2.14) [sin sN(s) ] ' -T (s) cos ^(s)-sin sH (s) f ' (s) - Ar (s) n ( s) =0 , 

/\ A A 

(2.15) [sin sH (s) ] ' +T (s) sin ^(s)+sin sN (s) ' (s) +Ar (s) v (s) =0 , 

A A A A 

(2.16) [sin sM ( s) ] ' - Z ( s) cos 4>(s)+sin s [v (s)H(s) -n (s)N(s) ] =0 . 

We can combine (2.14) and (2.15) to obtain 

A A 

(2.17) {sin sfsin \p (s) N (s) +cos ^ (s) H (s) ] +Ar (s) 2 /2 } '=0 , 
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A A 

(2.18) {sin s[-cos ^(s)N(s)+sin ^(s)H(s)]}' + 
+ T (s) +Ar (s) z ' (s) =0 


The material of the shell is homogeneous and non- 
linearly elastic if there are constitutive functions N, 

T, H, M, E such that 

/s 

(2.19) N ( s) = N(w(s)), etc., 

when (2. 3) -(2. 8) hold. 

We assume that these constitutive functions are thrice 
continuously differentiable on their common domain of 
definition 


(2.20) P = { w e R 5 : v>h | y+1 | , x > h|a+l| } 

where h is a given number in (0,1) , which may be inter- 
preted as half the thickness of the shell. The inequalities 
in (2.20) ensure that in a standard three-dimensional 
interpretation of our deformation variables the local ratio 
of deformed to reference volume never be zero. 

We require the constitutive functions to satisfy 
the monotonicity conditions : 


9(N,H,M) 
9( v, tv y) 


9 (T, Z) 

9( x, a) 


( 2 . 21 ) 


f 


are positive-definite. 
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( This is the strict form of the strong-ellipticity condition 
for axisymmetric deformations of axisymmetric shells. Cf . 
Antman (1978).) These conditions ens'ure,e.g. , that an 
increase in the bending strain y is accompanied by a 
corresponding increase in the bending couple M . We also 
use special forms of the following very natural inequalities 

(2.22) N + N >0 , T + T > 0 

V T TV 

We impose compatible growth conditions ensuring that extreme 
strains are accompanies by extreme resultants: 

(2.23) n (w) -*■ + 00 as v °° if x-h|l+o|has a positive 

lower bound and if ri, y,o are bounded, 
N ( W ) -*• - ® a s v -*■ h|l+y|if t is bounded above and 
if n,y,o, are bounded, 

T (w) -*• + °° as t ■+ 00 if v-h 1 1+y | has a positive 
lower bound and if n,y»o are bounded, 

T (w) -*■ — °° as v -*• h | l+o | if v is bounded above 
and if n»y»y are bounded, 

H (w) -*■ ± 00 as n ±°° if (v,x,y,o) lies in a 

compact subset of { (v,t ,y ,o) :v>h| l+o| , 
x>h | l+o | } 

M(w) -*■ t 00 as y -► ±°° if v >h | l+o | and if (T,ri,y) 
lies in a compact subset of {(T,n#y)j 
x>h | l+o | } 
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Z (w) -*■ + » as a -*■ ± » if x> h|l+a|and if (v,n,y) 
lies in a compact subset of ((v,n,y): 

v >h | 1+a | } . 

We finally require that the shell be isotropic and 
have a stress-free natural state : 

(2.24) n **• H (w) is odd, n h- N (w) , T (w) ,M (w) , Z (w) are even, 

M (w) = 0 = Z (w) if y = 0 = a , 

N a ( ~ )= % ( ~ )=T o ( ~ )=T y ( ~ )=0 if h = y = a = 0, 

N(v,x,0,y,a) = T(r,v,0,a,y) , 

M (v , r , 0 , y , a) = Z(T,v,0,a,y) , 

N (1,1, 0,0,0) = T (1,1, 0,0,0) = 0 . 

Our boundary conditions essentially require that the defor- 
mation be regular at the poles and that there be no net rigid 
displacement in the z-direction : 


(2 . 25a, b) 

r (0) = 0 = r (tt) 

(2 . 26a, b) 

n ( o ) = 0 = r» ( TT ) 

(2. 27a, b) 

1M0) = 0 , \p (n) 

(2.28) 

•TT 

z (s) ds = 0 

0 
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Substituting (2.25a) into the integral of (2.17) we obtain 

A A 

(2.29) sin s[N(s)sin i)j(s)+H(s)cos \p (s) ] +Xr (s) 2 /2 = 0. 


If N and H are bounded and if X > 0, 

then (2.27) implies that (2.25b) holds. If we further 

assume that t is bounded, then (2.29) and (2.27) imply 

A A ^ 

that H ( 0 ) = 0 = H ( 7T ) . The properties of H then ensure 
that (2.26) holds. By comparing the integrals of (2.18) 
over [0,s] and [s,tt] we readily obtain 


rrr 




(2.30) 


T ( s) ds = -X 


r (s) z ' (s) ds 


Note that the integral on the right side is the signed area 
between the curve (r,z) and the k-axis. 

Our boundary value problem consists of the strain 
configuration equations (2. 5) -(2. 8), the equilibrium 
equations (2.14)- (2.16) the constitutive equations (2.19) 
and the boundary conditions (2. 25) - (2 . 28) . A reaular 
solution of the remaining equations would automatically 
satisfy (2.25b) and (2.26). 

The shell theory we employ was originated by the Cosserats 
(1909). (Cf. Naghdi (1972) and Libai & Simmonds (1983).) Our 
formulation was foreshadowed by several papers of Reissner, 
typical of which is his paper of 1963. The variables we 
employ, and the very simple form of constitutive equations they 
induce, were introduced by Antman (1971) . 
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§3. TRIVIAL SOLUTIONS . 

We now seek solutions of our boundary value problem 
in which the shell remains spherical, unsheared, and uni- 
formly stretched, so that v = t = k(const), n = 0, ^(s) = s. 
Thus y = 0 = a . Under these conditions, (2.22) reduces 
(2.14)- (2.17) to 

(3.1) N(k,k, 0,0,0) = T(k,k, 0,0,0) = -Xk 2 /2 . 

Then (2.22) and (2.23) imply that k N(k,k, 0,0,0) + 
+Xk 2 /2 strictly increases from - °° to 00 as k increases 
from 0 to °° provided that \ ^ 0 . Thus (3.1) has a 
unique solution for k, denoted k(X) with 

(3.2) [0,°°) 3 X ^ k(X) e (0,1] 

thrice continuously differentiable and strictly decreasing . 
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§ 4. THE LINEARIZED EQUATIONS . 

We adopt the convention that if R is any constitutive 
function (such as N, N^, ), then 

(4.1) R° (X) = R(k(X) ,k(X) ,0,0,0) . 


The linearization of our boundary value problem about the 
trivial solution is equivalent to the following boundary 
value problem for the linearized variables v lf m , : 


f s 

(4.2) rj(s)= { v x ( t ) cos t- [m (t) +k^i (t) ] sin t} dt = 

o 


= xi (s) sin s , 


(4.3) N°[(V!sin s)'-TiCOS s] = [ (N® +H° +Xk)n i + 

+ (N°k+Xk 2 /2) i|>i] sin s , 

(4.4) H°(msin s) ' = - (N 0 +N°+Xk) ivi + Ti)sin s + 

n v t 

+(Xk 2 /2) (i|;isin s) ' , 


(4.5) M®[(i|>i'sin s)'-^i/sin s] + (M° -m£ Ji^sin s = 

= - [H°k+Xk 2 /2] msin s , 

(4.6) Tii(O) = 0 = m (it) , V»i(0) = 0 = Tpi (-rr) , 

■7T [IT 

(4.7) Zi(s)ds = {vi(t)sin t+ [m (t) +k^i (t) ] cos t}dt = 0 . 

Jo Jo 
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In deriving this system, we relied crucially on (2.24). 

From this system, we wish to extract a system for 
rii and ipi alone. To accomplish this efficiently we write 
(4.3), (4.4), and the derivative of (4.4) in the form 

(4.8) (visin s)'-Ticos s = f, (Vi+Ti)sin s = g, 

(Visin s) 'tvicos s= h , 

from which we immediately deduce that 

(4.9) g = (h-f)tan s 

We write out this equation, replacing (4>i'sin sj'-^x/sin s 
by it expression from (4.5) and supplement the resulting 
equation with (4.5), (4.6) to get the system 

(4.10a) (Lni) (s) + Arii (s) sin s = -aipi (s) sin s , 

(4.10b) (Lip x ) ( s) + -B^!(s)sin s = -brix(s)sin s , 

(4.11) Oi(0)= 0 = mU) , ip! (0) = 0 = (xr) , 
where 

(4.12) (Lu) (s) = [u'(s)sin s] ' - u(s)/sin s , 

(4.13) A = (H°) _1 [%Xk 2 b+(N°+N°+Xk) (N°) _1 (N®-N° + 

+Xk+H°)] , 
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a = (H° ) -1 [%Xk 2 B+ (N° +N° +Xk) (N° ) -1 (kN° - 
n v T v T 

-kN° +%Xk 2 ) ] , 

B = 1 - M°/M° , b = [H°k+Xk 2 /2] /M° 

We now convert (4.10), (4.11) to a problem for a 

single complex-valued function 

(4.14) 4> = am + Sty i 

where a and 6 are complex constants to be adjusted so 

that (4.10a,b) can be combined into a single complex 
equation. Thus we require that 


(4.15) (aA+8b)m + (aa+BB)ij;i = C<)> = C(ani+B^i) 

where C is a complex constant to be determined. It 
follows from (4.15) that 

(4.16) aA + 8b = aC , aa + 8B = BC 
so that 

(4.17) aa 2 + (B-A) aB - bB 2 = 0 . 
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Using (4.15) we reduce (4.10), (4.11) to 

(4.18) (L<j>) (s) + C<(>sin s = 0 , <f>(0) = 0 = <f> ( tt ) 

which has a nontrivial solution if and only if 

(4.19) C = n(n+l) , n = 1,2, 

in which case 

(4.20) <f>(s) = 4> n (s) = const, P* (cos s) 

where P 1 is the associated Legendre function of the first 
n 

kind of degree n and order 1 . When (4.19) holds, 
system (4.16) has a nontrivial solution if and only if 

(4.21) g(X;n) = [A-n (n+1) ] [B-n (n+1) ] - ab = 0 

Note that A, B, a, b each depneds on the eigenvalue parameter 
X . For a given integer n there can be none, one, or many 
real solutions X of (4.21) . Noting that N°(X)= -Xk(X) 2 /2, 
we could readily fashion conditions on the constitutive 
assumptions to ensure any of these possibilities. (Cf. 

Antman & Rosenfeld (1978) for a discussion of this issue.) 

We limit our attention to positive solutions X of (4.21). 
Negative solutions correspond to shear instabilities. (Cf. 
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Antman & Carbone (1977).) 

For X > 0, (2.21) implies that b > 0. Thus (4.16) 
has no nontrivial solutions ( a, 6) with a = 0. We can accordingly 
adopt the special normalization that a = 1 . In this case, 
the roots 6 of (4.17) are 


(4.22) 


_ 3-A t / (A-B) 2 +4ab _ B-A± [A+B-2n (n+1) ] 

± ~ 2b 2b 


the second equality coming from (4.21). Thus- 6 + are real 
when (4.21) holds. In this case, 8 + are distinct if and 
only if 

(4.23) A+B ^ 2n(n+l) 

When (4.21) and (4.23) hold, (4.14) and (4.20) reduce 
to a system of the form 

(4.24) m + 8 + ^i = D + P^(cos s) 

If we take D + real and not both 0 , then the distinctness 
of S + ensures that (4.24) can be solved uniquely for m 
and . These solutions are proportional to P*(cos s) 

To find D + , i.e., to find the constants of proportionality, 
we normalize by setting 

(4.25) ^i(s) = P^(cos s) , Hi (s) = EP*(cos s) 
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Then a comparison of (4.10) and (4.11) with (4.18) 
implies that 

(4.26) A + a/E = n (n+1 ) = B + bE 


whence 

(4.27) 


n (n+1 ) -B 
b 


a 

n (n+1 ) -A 


The compatibility of these two expressions for E is 
ensured by (4.21) . 

Now we study the special case that 

(4.28) A + B = 2n(n+l), n = 1,2, 

Then (4.22) implies that B is unique with 

(4 . 29a, b) 2b8 = B - A , bB 2 = a 

Consequently (4.20) yields 

(4.30) m + 6^i = DP*(cos s) 

Note that (4.29b) implies that there are no such solutions 


if a < 0 
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First we treat the even more special subcase that 

(4.31) A = B = n(n+l) 

Then 6 = 0 , m = DP^(cos s) . Moreover (4.16) implies 
that a = 0. Since b > 0 , the alternative theorem 
implies that (4.10b) has no solution unless D = 0 and 
hence n i = 0 . In this case, = const P^(cos s) . Note 

that if we make the eminently reasonable assumption that 
B < 2 , then (4.31) cannot hold. Likewise, (4.31) cannot 
hold if a # 0 

Now we assume that (4.28) holds, but that 

(4.32) A ± B * n (n+1) 

Hence 6 # 0. Using (4.30) we convert (4.10b) to 

(4.33) Lipi + n(n+l)4>! = -bDP^(cos s)sin s , 

since (4.28) and (4.29) imply that B - bB = n(n+l) . 
Since b # 0 , the alternative theorem implies that (4.33) 
has no solution unless D = 0 , in which case (4.30) yields 

(4.34) hi = 
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In summary, the eigenvalues of (4.10) , (4.11) 
are solutions X of (4.21). We limit our attention to 
positive' eigenvalues . If B < 2 , then the corresponding 
eigenfunctions are 

(4.35) (ni,^i) = ( [n (n+1) -B] /b, l)P*(cos s) 

Now the function P^(cos s) has exactly n+1 zeros on 
[0,tt] , including those at 0 and tt , each of which is 
simple. Thus every nontrivial solution (m»^i) of (4.10) , 
(4.11) is characterized by the fact that m and if)! have 
exactly the same zeros . 
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§5. GLOBAL BIFURCATION . 

To show that our equations meet the hypotheses of 
the bifurcation theorem and uniqueness theorem, which form 
the basis of our analysis, we must transform them into 
alternative forms that enable us to control the polar sin- 
gularity manifested by the presence of sin s as a coef- 
ficient of the highest derivatives of our differential 
equations. We substitute (2.19) into (2 . 14 ) - (2 . 17) , 
carry out the differentiation of the leading terms, use 
Cramer's rule (justified by (2.21)) to solve the resulting 
equations for v'sin s, n'sin s, y'sin s, and finally 
force this latest version of the system into a form suggested 
by the linear equations (4.3) — (4.5) . Using (4.12) we 
obtain 

(5.1) (v sin s) ' = t cos s + (N^) ^ (N°+H°+Xk) n sin s + 

+ (N^) 1 (N°k+Ak 2 / 2 ) ( s ) sin s + n , 

(5.2) (n sin s) ' = - (H° ) _1 (N°+N°+Ak) ( v+t) sin s + 

+ (H®)‘ 1 (Xk 2 /2) [ (i(»-s)sin s] ' + h , 

(5.3) L ( ip-s ) = -B(^-s)sin s - b sin s + m , 
where 

(5.4) n = A _1 D 1 + ( v— x ) cos s- (N° ) _1 (N°+H° +Ak) n sin s - 

- (N°) _1 (N°k+Ak 2 /2) (^-s) sin s , 
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(5.5) 

(5.6) 

(5.7a) 

(5.7b) 

(5.8a) 
(5.8b) 
(5.8c) 
From (2. 

(5.9) 

(5.10) 


h = A -1 D 2 +ticos s+ (H°) _1 (N° +N° + Xk )(v+r)sin s 
- (H° ) -1 (Xk 2, /2) [ (\p— s) sin s] ' , 


sm s 


m = A D 3 + (^-s) 'cos s - - f~ s g + B ( ip- s)sin s 
+ bnsin s , 


A = det 

3(v,n»v) » 


is the determinant obtained from A by 


replacing its jth column with j 8 

Y ' 


a >, 


a = - (N^t ' +N^a 1 ) sin s+Tcos ijj-Ncos s+Hty'sin s+Xrri , 


8 = - (H^t ' +H a a ' ) sin s-Tsin ;|;-Hcos s-N^'sin s-Xrv , 


y = - (M^x ' +M a a ' ) sin s+Ecos i|;-Mcos s-(vH-nN)sin s 


5) -(2. 8) we find that 


Lr = (vsin s) 'cos \p+ (vsin s) (cos ip) ' - 

-(nsin s) 'sin ip- (nsin s) (sin ip) ' - t , 

(z'sin s) ' = (vsin s)'sin ij;+(vsin s) (sin \p) ' + 

+ (nsin s) 'cos i^+ (nsin s) (cos \p) ' . 
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We now introduce new variables (u!, U 2 , U 3 ) = u by 
(5.11a) ui (s) (sin s)^ = (Lr) (s) , 

(5.11b) u 2 (s) (sin s )* 5 = (z'sin s) ' , 

(5.11c) u 3 (s)(sin s )* 5 = L(^-s) (s) 


Let us assume that u is continuous. Then we can solve 
(5.11a,c), (2.25), (2.27) in terms of ui and u 3 . Now 

from (5.11b) we see that u 2 must satisfy 


(5.12) 


r 7T 

u 2 ( s ) (sin 

0 


s) ^ ds 


0 . 


We thus find that 


(5.13) 


z ' (s) = (sin s) 


- 1 


u 2 (t) (sin t) ^ dt 


= -(sin s) 


_ i 


r rr 


u 2 (t) (sin t) ^ dt 


This representation shqws that z satisfies the boundary 
conditions z'( 0 ) = 0 = z' ( tt ) , which are consequences of 
(2. 5) -(2. 8 ), (2.26), (2.27). We integrate one of the 

equations of (5.13) and then use (2.28) to evaluate the 
constant of integration. We obtain 


nr 


(5.14a) r(s) = 


K(s,t)ux (t) (sin t)^ dt = (Gux) (s) , 
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(5.14b) z(s) 


f tTr~ 1 y-l+x [0>sl (Y)1 

o sin y 

= ( Zu 2 ) (s) , 


y , l 

u 2 (t) (sin t)^dt dy 

' o 


(5.14c) iMs)-s 


(TT 

J 0 


K ( s , t) u 3 (t) (sin t) 



(Gu 3 ) (s) 


where Xrn e i is t ^ ie characteristic function for [0,s] 
and K is Green's function for L subject to Dirichlet data 


(5.15) 


, 1+cos s 
^ 1+cos t 


sin t 
sin s 


K(s,t) 




1+cos t sin s 
1+cos s sin t 


for t < s , 


for t > s 


(We could express (5.14b) in terms of a Green's function, 
but doing so would interfere with our analysis. ) We also 
record that 


(5.16a) 

r' (s) = 

K (s,t)Ui (t) (sin t) *dt = (GUi) ' (s) , 



0 

(5.16b) 

(s) -1 

k 

= K g (s,t)u 3 (t) (sin tpdt = (Gu 3 ) ' (s) 



0 


where. 
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(5.17) 


K s (s,t) = 


z' is given by (5.13) . From (2.7), (2.8) we then obtain 

(5.18a) v = (Gu x ) ' cos (s+Gu 3 ) + (Zu 2 )'sin (s+Gu 3 ) , 

(5.18b) n = “ (Gux) 1 sin (s+Gu 3 ) + (Zu 2 ) ' cos (s+Gu 3 ) , 

(5.18c,d) t = (GUi)/sin s , a = [sin (s+Gu 3 ) /sin s] - 1 

We now convert our differential equations (5.1)- 
(5.3) into integral equations for u : We replace the 

left sides of (5.9), (5.10), and (5.3) with (5.14) . We 
next substitute (5.1), (5.2) into the right sides of (5.9), 

(5.10) . Then we replace the variables appearing on the right 
sides of (5.9), (5.10), (5.3) with their representations 
from (5 . 14) - (5 . 18) . Finally, as a precautionary measure, 

we subtract from the right side of the resulting form of 

(5.11) its mean value, a move suggested by (5.12). We 
denote the resulting system of integral equations by 


1+cos s 
1+cos t 


sin t 

7 — ~ 


sin 


for 


t < s 


-h 


1+cos t 1 

1+sos s sin t 


for 


t > s 


(5.19) 


u = f (A,u) 
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Let 

(5.20a) U = { u £ C 0 ( [ 0 , tt] ) 3 : (5.12) holds.) 

We wish to show that f(X,‘) is a compact mapping from 

(5.20b) U fi{ u : w e V } = E 

into itself ( V is defined in (2.20)). 

Let w -*■ R(w) be any constitutive function. If 
R e C 1 ( V ) , then 

(5.21a) R (y ) = R(x) + (x,y) (y^-x^ 

with 

f 1 

(5.21b) R i (x,y) = R^ (ty+(l-t)x) dt 

o i 

Here and below, twice repeated indices are summed from 1 to 
5. If R e C 2 ( V ) , then 

(5.21c) R(y)=R(x) + R (x) (y . -x . ) +%R. . (x,y) (y . -x . ) (y -x . ) 

~ ~ x -*■ - L J ~ ~ ■*- 1 3 J 

with 

(5 . 21d) R ii < x 'y) = R w w dt 

J 0 i 3 
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We let R (w) and R (w) be any values of R. and R. . 

w. ~ w.w. ~ L 1 in 

ii] J 

in which the arguments consist solely of entries from w. 

E.g., R (w) could be taken to be R.((T,v,0,a,u), 

W • ~ 1 

1 

(v,T,n,y,a)) . 

Since we whall use only the continuity of { R ,R } 

w i* W i W j 

and not the specific functional forms, we allow their 
meaning to vary even in the same expression. 

In the spirit of (5.21) we define the continuously 
differentiable functions C and S by 

(5.22a) 1 - cos 0= 0 2 C(0) , 

(5.22b) 0 - sin 0 = 0 3 S(0) . 

Then 

(5.23a) cos \p- cos s = - (ij;-s)sin s + U , 

(5.23b) U = - (i^-s) 2 C (i^-s) cos s + (\p-s) 3 S ( ^»— s ) sin s f 

(5,23c) sin iJj - sin s = (^-s)cos s + V , 

(5.23d) V = - ( s ) 2 C (tj'-s) sin s - (\^-s) 3 S ( ip— s ) cos s 


Thus 
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(5.23e) t ' sin s = (v-t)cos s - (^-s) (vsin s+ycos s) - 

- nsin s + vU - nV , 

(5.23f) a'sin s = (y-a)cos s + ip 1 [U- ( ip— s ) sin s] 

Let 

(5.24) w* = (x,v,0,a,y) 

By using the formulas (5 . 21 ) - (5 . 23) we can convert a , 
given by (5.8a), to 

(5.25) a = T (w) [U- (^-s) sin s] + % T y 2 cos s - 

-[N (w*) { v— t ) +N (w*) (y-c) +* 5 ^ . (w* ,w) (w i ~w i ) 

* * ★ 

(w.-w )+N (v-T ) (w. -w. ) +N (y-a) (w. -w )] cos s 

J J TW ^ 1 X UW^ X X 

+ (w) [ (^-s) (vsin s+ycos s)+ysin s-vU+yV] - 

-N^ (w) \p ' [U- (^-s) sin s] + H^yip’sin s + Ary 

We obtain similar representations for 8 and y . Note 
that a + !N v (w*)(v-T) + N u (w*)( l i-a)]-n(H I1 *'sin s + Xr) consists 
of sums of product of continuous functions of (w,i^-s) with 
quadratic functions of v -t, y, y-a , ip-s . By a lengthy 
computation we find that the terms in f consist of products 

of continuous functions of (w,^-s)with quadratic functions 

-i' _y — y 

of (sin s) 4 (v-x) , (sin s) 4 y , (sin s) 4 (y-a) , and 
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with a number of other such functions of (w, ^-s) . Let 
us examine (sin s) 4 (y-a), which (5.23c), (5.14c), (5.16b) 
reduce to 


(5.26a) 


(sin s) 4 [i|>'(s) - l - (ij; ( s ) -s) cot s - V(s)csc s ] 

,7r - 5 / 

J(s,t)u 3 (t) dt - (sin s)" /4 V(s) 

o 


where 


(5.26b) 


J (s, t) 


r 


H 


(1+cos s) 2 f sin s ‘ 
1+cos t ^sin t 



(sin t) 


/ 4 


t < s, 


< 




-h 


(1-cos 

1-cos 



( . . •N 

sm t 
>sin s. 


9 /4 


(sin t) 


/ 4 t > S 


The Arzelci-Ascoli Theorem supports a straightforward proof 

that the kernel J generates a compact and continuous 
operator from C°([0 ,tt]) to itself. Note that the kernel 
(sin s) _3s K ( See (5.17)) does not. It is the cancellation 

o 

of terms in J that is responsible for this compactness. 

Now (5.23d) implies that 


_ 5 / _ 1 / 

(5.26c) -(sin s) 4 V = [(sin s) 8 (^-s) ] 2 C (^-s) + 


_ 5 / 

+ [(sin s) ' 1 2 (^-s) ] 3 S (i^-s) cos s 


Since the Arzela-Ascoli Theorem implies that the integral 

_ 5 , 

representation for (sin s) /12 (^-s) obtained from (5.14c) 
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in a compact mapping of u , it follows that (5.26c) 
defines a compact mapping. Since compositions of continuous 
and compact mappings are compact, we find that (5.26c), 

_ i / 

(5.26a), and (sin s) /2 (y-cr) 2 correspond to compact 
operators taking C° ( [0,tt] ) into itself. We emphasize that 
the appearance of terms such as (y-a) in f is a consequence 
of the isotropy conditions. 

By extensive use of arguments such as these. We 
find that f is a compact and continuous mapping of E 
into itself. Moveover, the form of f shows that it is 



let there be an odd positive integer 


such that 
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. a j 

(5.29) ? g(X .u) = 0, j=0 , . .k-1 , -2— g (X ,u) + 0 . 

9X J n ~ ax K n ~ 

Then bifurcating from the point (X n ,u 0 (X n )) on the trivial 

+ 

branch of solutions are two connected sets C~(X n ) of non- 
trivial solution pairs (X,u) of (5.19) that enjoy at 
least one of the following properties : 

(i) C - (X n ) cannot be confined to any compact subset 

of E 

(ii) The closure of C ± (X n ) contains another point 
( c,uo (?) ) where z, satisfies (4.21) 

Suppose, moreover, that X n is simple , i.e. , that k=l 

in (5.29). Let v n be a normalized eigenfunction of (5.27). 

Let P (X ) be the projection of U onto span(v_). Then there is 
— n k ~n 

a number T > 0 and there are continuous functions < : 

n 

[-X, X] h- r, : [-X,X] h- [i-p (X R ) ] u with < n (0) = 0 , 

y n (0) = 0 such that (X n +K n (e), u 0 (X n +< n (e)) + e[v n +y n (e)]) 

is a solution pair of (5.19) for | e | < X . Moreover, 
there is a neighborhood W n of (X n ,uo(X n )) such that if 

(X,u) is a solution pair of (5.19) lying in > then 

either (X,u) = (X ,u 0 (X)) or else there is a n e e [-X,X] 

such that (X,u) = (X +< (e) , u 0 (X_+<_(e)) + e[v +y (e)]\ . 
n n ~ n n ~n ~n 

Note also that (5-14) — (5-18) imply that w cor- 
responding to any solution of (5.19) is continuous and in 
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particular \p e C 1 ([0,t t] ) . Moreover 

(5.30) lim [v(s)-t(s)] = 0 , lim [y(s)-a(s)] = 0 

S^-O/TT •S->0,7T 

It follows from (5.1) -(5. 3) that u generates a classical 
solution to our boundary vlaue problem. In fact, the re- 
quirement that the constitutive functions N, T, H, M, Z 
be thrice continuously differentiable on p implies that 
the strains w generated by u are twice continuously 
differentiable on the open interval ( 0 , tt ) . To determine 

the boundary behavior of n’ , which is needed in the 
sequel, we divide (2.29) by sin 2 s and use (2.19) and 
(5.21) to obtain 

(5.31) N(w(s) ) S gY n ^p ~ +H^ (w) s ^ S s cos ’Ms) +Xt(s) 2 /2 = 0 

on ( 0 , tt ) . Since r\ satisfies (2.26), we can let s ■+• 0,ir 
in (5.31) to conclude that n is differentiable on [ 0 , it ] 
To show that n e C 1 ([0 ,tt]) we use the expression for n' 
obtained from (5.2) and examine its limits as s -*• 0,it 
This process is simplified by the observation that the 
limiting forms for n 1 can equivalently be obtained by 
substituting (2.19) into (2.15), carrying out the 
differentiations, dividing the resulting equation by sin s 
and exploiting (2.24) and (5.30) to show that the limits 
n ' (0) and n’('n') agree with those found from (5.31). 
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essentially equivalent version of Theorem 5.28 can 
be obtained by usmq the theory of Fitzpatrick & Pejsachowicz 
(1985), which generalizes that of Crandall and 3abinowitz. The 
former theory does not require the use of Green functions. 

The role of isotropy still remains paramount. 
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§6 . NODAL PROPERTIES OF SOLUTIONS 

The determination of nodal properties of 
solutions of our boundary value problem is a more subtle 
question than that for the boundary value problem for 
plates. The main reason for this difference is that in the 
plate theory (of. Antman (1978)) n is the product of a 
certain function with sin i p and so the nodal properties of 
n are determined directly from those of ij/ . On the other 
hand, in our shell theory v , x , n , ^ are all 
coupled by nonlinear differential equations, which lack an 
obvious nodal pattern. To approach the question for shells 
we first need a candidate for a nodal pattern that is 
preserved on solution branches. 

The solutions of the linearized problem, treated in 
Section 4 , have the property that m and ipi have the 

same zeros, which are simple. The perturbation solutions 
of Section 8 also have this feature. Thus we are motivated 
to study whether the number of simultaneous zeros of n and 
ip - s is fixed on each bifurcating branch. 

We immediately confront a serious difficulty. Let 

^ n ' n = 1, 2, ••• be the collection of all real-valued 
functions p in C 1 ( [0 ,tt) ] having zeros at 0, ir and 
having exactly n - 1 zeros on (0,ir), each zero being 
simple. Then 5^ is open in C 1 ([0 ,tt)]. Given any 

function in , it follows that all sufficiently nearby 

functions enjoy the same nodal properties, i.e., have the 
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same number n + 1 of zeros, each of which is simple. 

Moreover, the boundary 95^ consists of functions 

having at least one double zero on [ 0 , tt ] . These observations 

support the nodal theory introduced by Crandall & Rabinowitz 

(1970) for bifurcation problems governed by nonlinear 

Sturm-Liouville equations. To study pairs of functions with 

simultaneous zeros, we let T be the collection of all 

n 

complex-valued functions in [C 1 ( [ 0 , tt] ) ] 2 having zeros at 

0, tt and having exactly n - 1 zeros on ( 0 , m ) , each zero 

being simple. But is not open in [C 1 ( [0 ,it] ) ] 2 and 

there is no useful characterization of 9 T . To 

n 

circumvent this difficulty, we must exploit the fact that 
(rif 4>~s) is not merely an arbitrary element of T , but 
is also a solution of a boundary value problem. 

Our strategy is as follows. We first describe an 
auxiliary condition ensuring that if there is a solution of 
our boundary value problem for which n , r \ 1 , \p -s, 

(iJj-s) ' all vanish at a point So in [ 0 , rr ] , then v-t 
also vanishes at So • Next we prove a basic uniqueness 
theorem stating that if there is a solution of our boundary 
value problem to which v-t, ri , ip - s , (ij/-s) all 

vanish at a point So in [ 0 , tt ] , then the solution is 
trivial. This result does not enable us to assert that on a 
nontrivial solution branch (r\, ip-s) remains in T n if 
(n, ip— s ) is in at any solution pair on the branch. 

To remedy this difficulty, we set 
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(6.1) n(s) = p(s) cos 0(s), 

i^(s) - s = p(s)sin 0 (s) . 

If n(s 0 ) = 0 = ^(s 0 ) - so / then p(so) = 0. Thus the zeros 

of the real-valued function p are the simultaneous zeros 

of ri and ip - s . If n and ip are part of a nontrivial 

solution of our boundary value problem for which the auxiliary 

condition holds, we use the uniqueness theorem to prove that 

p can be defined so that it belongs to C 1 ([0 ,tt]). We 

can accordingly determine where p e S . In particular, 

if at a point on a nontrivial solution branch pc S , 

then p remains in S n until it has a double zero ( by 

virtue of our discussion of 95 ) . 

n 

By the remarks at the beginning of the last paragraph, 

we conclude that the solution is trivial if p has a double 

zero ( and hence ri and \p - s have simultaneous double 

zeros ) provided the auxiliary condition in met. Let X n be 

a simple solution of (4.21) and let C ± ( X ) be the two 

n 

nontrivial branches bifurcating from (X n ,uo(X n )). We wish 

to show that on C " ( X ) the function p remains in 

n 

as long as the auxiliary condition is met. To do this, 

we must show that p admits a linearization and enjoys 

properties like those stated in the last part of Theorem 5.28. 

We thereby demonstrate that on C ± ( X n ) , the function p 

inherits its nodal properties from those of its linearization 

and preserves them as long as the auxiliary condition is met, 

i.e., on C~ ( X ) , p e S as long as the auxiliary 

n n 
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condition is met. Thus on C ± ( X ) the function n and 
ip - s have n - 1 simultaneous zeros on ( 0 , tt ) as long as 
the auxiliary condition is met. 

We now proceed to carry out this program. We assume 
throughout the ensuing discussion that (X,^,w) generates a 
solution pair for our boundary value problem. 

Suppose that there is an s 0 £ ( 0 , tt ) such that 
n(s 0 ) = ri'(so) = iMs 0 ) - so = ^’(s 0 ) - 1 =y(s 0 ) = 0 . Then 

(2.24) ensures that (2.15), (2.19), (2.29) reduce to 

(6.2) N + T + Xx(s 0 )v(s 0 ) = 0, 

(6.3) N + h Xx (s 0 ) 2 = 0 

where the arguments of N and T are (v(so) , x(s 0 ), 0, 0, 0). 
We know from Section 3 that (6. 2), (6. 3) admit a solution 
of the form v(s 0 ) =* x(s G ) = k(X) and that this is the only 
solution of these equations for which v(so) =x(s 0 ). By 
virtue of (2.20) the variables v , x , X are confined 
to (h,°°) x (h, 00 ) x (0,°°) . Let G be the subset of this 
region for which the only solution of (6.2), (§.3) is 
v(s 0 ) =x(s 0 ) = k(X). We then have the following result, 
which is virtually a tautology. 

6.4. LEMMA . If n(s 0 ) = n'(So) = iMs 0 ) - So = f (s 0 ) - 1 
= 0 and if (v(so), x(sj) ,X) e G , then 
v(s 0 ) = x(s 0 ) = k (X)^ 
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In the next section we shall study the nature of G 
the "auxiliary condition", alluded to above is just the 
requirement that ( v ( s 0 ) , x (s 0 ) , X) e G . Our basic 
uniqueness theorem is 

6.5. Theorem. If there is an s 0 e [0 ,tt] such that 
(6.6) v(s 0 ) ~ k(A) = x(so) - k(A) = n(s 0 ) = <Ms 0 ) - s 0 = 

= $ ' ( s o ) ~1 = 0 , 

then the s'olution is trivial . 

Proof. Let 

(6.7a) v = (v : ,v 2 ,v 3 ,v 4 ,v s ) 

= ((v-r)sin s, (v+t-2k) sin s, nsin s, ysin s, \p- s) , 

(6.7b) V = j Vi | + |v 2 | + [ v 3 j + | v «+ | + | v 5 | 

Then we can combine (5.1) — (5.3) # (2.8), and (5.23) to 
obtain the system 

• -1 
(6.8a) v x = v(v 5 sin s-U) 4v 3 + n (v 5 cos s +V) + A D x , 

I 

(6.8b) v 2 = (v x +v 2 )cot s - v(v 5 sin s-U) - v 3 - 


/ 


-n(v 5 cos s+v) + a ‘*’D 1 
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(6.8c) 


v 3 


v 3 cot s + 



f 


(6 . 8d) 


■ 

V 4 


V 4 COt 


S + 



D 3 


/ 


' "I 

(6 . 8e) v 5 = v 4 (sin s) 

If so e (0 ,tt) and if v satisfies (6.6) and (6.8) 
then v = 0 , because the right-hand side of (6.8) is a 
well-behaved function of s and v . Consequently the 
initial-value problem has a unique solution, which is the 
trivial solution v = 0 

Thus we need only concern ourselves with the cases 
So = 0, it, in which the right-hand side of (6.8) is 

singular in s . We just treat the former case. Throughout 
our discussion we let v represent an arbitrary, given 
solution of our boundary value problem ( known to exist by 
the results of Section 5 ) which we subject to the 
additional requirement that it satisfy (6.6) . In the 
analysis to follow, we let c denote a typical positive 
constant depending on the fixed but arbitrary solution v . 

The meaning of c can change from place to place. We tacitly 
restrict s to lie in an interval of the form [0,e] where 
z is small and positive. 

Since sin s ^s cos s on [0,e] , we find that 


r s 


(6.9) 


-h h h 

(sin t) ^dt < 2(sin s) vcos s < c(sin s) 
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Using this result in (5 . 13) - (5 . 18) and using the fact that 
| u | £ c , we obtain after careful analysis that 

(6.10) | v — t | + | n | + |y-a| £ c(sin s) 1 * , 

|ijr-s| £ csin s , |y| + |cj| £ c 

Since (5.23) implies that 


( 6 . 11 ) 


a| £c|v 5 |(sin s) 


-1 


we find from (5.25) and related expressions that 


( 6 . 12 ) 


a| + | 6 | + | y 1 < c(sin s) 


Let us now examine (6.8a) . From (6.10) we 

immediately obtain 


1 — 1 — 

(6.13) |vi-A Di| £ c(|v 3 |+|v 5 |) £ cv £ cv(sin s) ^ 

Note that (2.21) and (2.23) imply that 

(6.14) A -1 £ c • 

Accounting for (2.24) by using the notation introduced after 
(5.21) we obtain 
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(6.15) 


Di = aH M - aH M n 2 - (BM -yH )N n + 

n y yn nn y ' y nn 


+ BN M n ~ yH (N n+N y+N a) , 
y nn ' n yn yy yc ' 


(6.16) 


A -1 - (N H M ) -1 =Q=Qn+Qy+Qa 
v n y n y a 


Inequalities (6 . 10) - (6 . 12) , (6 . 14) - (6 . 16) imply that 


(6.17) 


A ^Di-atN^) *| £ cv(sin s) ^ 


If in (5.25) we replace N (w*) with N (w ) +N (w.-w.) 

r v ~ v - vw/ 1 i 7 

and replace N (w*) with N n+N y+N a , then we obtain 

k ^ y y y ^ 

from (6.10), (6.11) that 


(6.18) 


a +N (v-t)cos s < cv(sin s) 


->5 


Inequalities (6.13), (6.17), (6.18) imply that 


(6.19) 


1 — Jj 

; Vx+Vicot s| < cv(sin s) 


A virtually identical analysis yields 


( 6 . 20 ) 


V 2 — V 2 cot s £ cv(sin s) 


-h 


The treatment of (6.8c) 'is just slightly different. 
Using the same ideas as those leading to (6.17) we obtain 
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(6.21) | A 2 | £ c | 6 | + cv(sin s) ** 

Since = H ti, = H^n by (2.24), equation (5.8) 

yields 

(6.22) | 6 1 £ cv(sin s) ** + | T+N+Avx | sin s 
Equation (3.1) implies that 

(6.23) T + N + Xtv = X = X (v-k) + X (x-k) 

V T 

= j 5[X v (vi+v 2 )+X v (-vi+v 2 ) ] (sin s) 1 
Combining (6.8c) , (6 . 21) - (6 . 23) , we obtain 

(6.24) |v 3 -v 3 cot s| < cv(sin s) ^ 

We remark that the only place that the variable v 2 appears 
in (6.8a,c,d,e) is in the right side of (6.8c). Similar 
techniques produce 

(6.25) |v 4 - vscot s| <_ cv(sin s) ^ 

Combining (6.19), (6.20), (6.24), (6.25), (6.8e) we 

obtain 


v' < 



£ vcot s+cv(sin 



/ 


(6.26) 
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the first inequality being a standard consequence of the 
absolute continuity of v . Since (6.6) implies that 

(6.27) v(s)/s -*■ 0 as s -*• 0 , 

we can invoke the standard uniqueness theorems associated 
with differential inequalities (cf. Hartman (1964 , Sec. III. 6)) 
to conclude that v = 0 . g 

It is important to note that if we were to insert a 
constant exceeding 1 in front of vcot s in (6.26), then 
we could not conclude that v = 0 . That the coefficient 

in front of vcot s is 1 is a consequence of the isotropy 


condition. 
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We are now ready to study the polar coordinates p 
and 0 introduced in (6.1). We state their properties in 
a series of lemmas. 

6.28. Lemma . Let the boundary value problem have a non- 
trivial solution pair with (v,x,X) taking values only in 
G and with p and 0 defined in any manner consistent 
with (6.1). Then the zeros of p on [0,ir] are isolated . 

Proof. If not, they would accumulate at a point s 0 £ [0 ,tt] , 
at which *n and ip -s would consequently have simul- 
taneous double zeros. By Lemma 6.4 and Theorem 6.5, the 
solution would then be trivial. 0 
Now (6.1) implies that 

(6.29) tan 0= (^-s)/n . 

Suppose that the solution is not trivial and that (v,x,X) take 
values only in 6 . Let Si and s 2 be a pair of adjacent 

zeros of p If n vanishes nowhere on (si,s 2 )/ then 

(6.29 ) implies that tan 0 is continuously differentiable 
on (si, s 2 ) and we can define 0 itself (mod it) to be 
continuously differentiable on (si, s 2 ). If n vanishes at 
S 3 £ (si , s 2 ) , 


then we define 
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(6.30) 0{s) 


tt/2 + 


' s fn(t) ir (t)-l)- t »(t)-11V(t)} dt|rood ir) 
^S 3 T) (t) 2 + (t) -t] 2 


for s e ( s i , s 2 ) • If S 4 is a zero of p and if n' (s 4 ) 

/v 

and ^'(s 4 ) - 1 do not both vanish, then I'Hopital's rule 
enables us to deduce from (6.29 ) that 

(6.31) tan 0(s) ^ '^sl) ^ as s s 4 , 

the limit existing as an extended real number. (If n'(s 4 )= 

= 0 = ip' (s )-l, then we could use Lemma 6.4 and Theorem 6.5 
to show that the solution is trivial when the condition bn 
G holds.) Then 0(s 4 ) is defined mod tt . By selecting 
the value of 0 at such zeros of p appropriately, we can 
define 

(6.32) p = ncos 0 + (^-s) sin 0 
Hence Lemma 6 . 4 and Theorem 6 . 5 imply 

6.33. Lemma . Let the boundary value problem have a non- 
trivial solution pair with (v,n»X) taking values only in 
6 • Then p and 0 satisfying (6.1 ) can be chosen to 

be continuous on [ 0 ,tr]. 


Since ti,^ e C 1 ( [0,tt] ) , the constructions we have 
made above ensure that p and 0 are continuously differen- 
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tiable where p does not vanish provided the hypotheses of 
Lemma 6.33 hold. Differentiating (6.1 ) we obtain 


(6.34a,b) p' = n’cos 0 + (ty'-l)sin 0 , 

p0' = -n'sin 0+ (ip'-l)cos 0 

where p does not vanish. The continuity of the right 
sides of (6.34) then yields 

6.35. Lemma. If the hypotheses o.f Lemma 6.33 hold , then 
P' , P<3' e C° ( [ 0 , tt] ) . 


We require the following refinement of this lemma : 


6.36. Lemma. If s 0 is a zero of p and if n'(s 0 ) 

and ty'(s 0 )-l are not both zero , then p(s D )0 , (s o ) = 0. 

Proof. The Mean Value Theorem yields 


(6.37) n (s) = 


f z 1 


' o 


n ' (as+ (1-a) s 0 ) da 


(s-s 0 ) , 


r 


iMs) - s = 


V (as+ (1-a) so ) da 


(s-s 0 ) 


On any small one-sided neighborhood of s 0 / (6.34b) thus 

yield 
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(6.38 ) 


±p (s) 0 ' (s) 


y (s) 


ri ' da-n ' (s) 


yda 


[ 

f 

n ' da] 2 + [ 

\ 

yda] 2 

1 * 

o J 

o 


where the arguments of n ' 
as + (l-a)so . Since ri ' 
on [0,7 t] , we find that 


and y in the integrands are 
and y are uniformly continuous 
[ y ' (as+ (1-a) So ) da -*■ n'(s 0 ) as 


J o 

s -*■ so , etc. Thus in the limit as s + So , the numerator 
on the right side of (6.38) approaches 0 , while the 
denominator approaches a positive number . 0 


In view of Lemma 6.35 we can examine whether 
p e S . By the introductory remarks of this section, func- 
tions p in 9S have a double zero. But this fact alone 

n 

does not guarantee that n and \p-s have a simultaneous 

double zeros, as an examination of (6.1) and (6.34a) shows. 
But Lemma 6.36 coupled with (6.34) implies that at a 
double zero of p the functions n and \p-s have a simul- 
taneous double zero, a result that enables us to invoke 
Lemma 6.4 and Theorem 6.5 . Marshalling these facts in 

the manner proposed in the introductory remarks in this 
section, we obtain 


6 3 9. Theorem. Let K be a connected set of nontrivial 
solution pairs with the property that (v,t,A) takes values 
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only in q for each solution pair in k . If for one 
solution pair in K the function p e s n / then p 
remains in S n for every solution pair in k . Thus n 
and ip-s have exactly n+1 simultaneous zeros on k 

Let J n be a simple root of (4.21 ) (so that (5.29 ) 
holds with k = 1 ) . Then the last part of Theorem 5.28 
together with the representations (5.14c ) and (5.18b) 
imply that on the branch C( X n ) near the bifurcation point 
( X n , 0 ) , X , n , ip have representations of the form 

(6.40) X - X n + < (c) , n = em + o(e), 

- s = cilii + o(e) 

as e 0, where < is continuous and c(0) =0. If we 
substitute (6.40 ) into (6.29 ) and use (4.25 ) we 
obtain 

(6.41 ) tan 0 = /hi + o(l) = E _1 + o(l). 

Thus we may take 

(6.42 ) G = arctan (E _1 ) + o(l). 

Then (6. 32) implies that 

(6.43 ) P = e Pl + o (e) = e [ Pl +o (e) e" 1 ] 
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with 

(6.44 ) Pl = (1+E ) _l5 [Eru +tyi] 

= (1+E ) ** P 1 (cos s) . 
n 

It follow that Pi £ S and that p e S for e 

n n 

sufficiently small. ( Incidentally, we show in the next 
section that (v,x,A) generically belongs to G for e 
sufficiently small essentially because v - t= 0 on the 
trivial branch. ) We accordingly have 

6.45. Theorem. Let 7 be a simple root of (4.21). 

n 

On the maximal connected subset of the closure of q~ (X ) 

containing ( uo ( T n ) ) for which (v,t,A) takes values 
only in G , ri and -s have exactly n + 1 simul- 

taneous zeros. 
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§7. PROPERTIES OF G. 

We now study (6.2), (6.3), which we rewrite as 

(7.1) N(v,x) = -^Xx 2 , 

(7.2) N(v,x) + T(v,t) = -Xvx , 

V 

from which we obtain 

(7.3) T(v,t) = Xt (*st-v) 

Here we have suppressed the last three arguments (0,0,0) 
of the functions N and T . The set G consists of 
all ( v , t , X ) in (h,») x (h,°°) x (0,°°) for which the only 
solution of (7.1), (7.2) is v = x = k(X) . We make a 

few remarks about (7.1), (7.2) in general before studying 

a specific example. 

Since > 0 by (2.21) and since N satisfies 

the growth conditions of (2.23), equation (7.1) can be 
solved uniquely for v in terms of t and X : 

(7.4) v = f ( t; X) , with 1 = f (1 ; 0) . 

Since f = - (N +Xt)/N ,we find that f < 0 if N > 0 

T T V T T 

and that f ( 1 ; 0 ) > - 1 if N (1,1) < N (1,1) 

T TV 

The isotropy condition (2.24) ensures that the 
curve in the (v,x) -plane defined by (7.2) is symmetric 
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about the line v = t . If T + N >0, if T + N 

v v ' 

as v 0 for fixed x , and if liminf T + N ^ 0 as 
v ->• oo for fixed t , then (7.2) can be uniquely solved 
for v in terms of t and A : 


(7.5) 


v = g(x;A) with 1 = g ( 1 ; 0 ) 


The results of Section 3 ensure that the graphs of f(*;A) 
and g(*;A) intersect on the line v = t . Since 


(7.6) 


= - 


(Ag+N T +T T ) 

vw 


it follows that g^ < 0 under our assumptions and, not 
surprisingly, that f^fx;!) = -1 when x = g(x;A) . Thus 
the graphs of f(*;A) and g ( • ; X ) cross transversally 

on the line v = x when f ^ -1 here, i.e., when 

x 


(7.7a) 


N v (x»x) - N t (x,x) - Ax ± 0 


Using (7.1) we can rewrite (7.7a) as a condition involv- 
ing only the constitutive function N : 


(7.7b) N v (x,x) - N t (x,x) + 2N(x,x)/x ^ 0 . 


Wherever (7.7) holds, there is a neighborhood about 
(x,x,A) lying in G . ( See the parenthetical remark 
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preceding Theorem 6.45.) 

There are several ways to analyze (7.1), (7.2) 

further by making general restrictions on N and T 
leading to suitable estimates. It is also illuminating to 
view (7.1) and (7.2) as each describing the inter- 
section of two graphs over the (v,x) -plane. Alternatively 
we could assume that N and T are derivable from a 
stored energy function { as in (7.9) below ), characterize 
the equations as Euler-Lagrange equations for a constrained 
problem, and study the extremization of a real valued 
function on a constraint curve. There are several ways to 
choose the constraint. It is more illuminating, however, 
to analyze a specific class of materials. 

Let A, B, C, a be positive numbers and let 

(7.8) W(v,t) = Act a t a + j sB(v+t) 2 + j sC(v 2 +t 2 ) , 

(7.9) N ( v , T ) = W v (v,T) - W v (l,l) , 

T(v,t) = W t (v,t) - W T (1,1) 

These equations give a reasonable model for material behavior 
in compression for h = 0 . When (7.8), (7.9) hold, 

equations (7.1), (7.2) reduce to 

(7.10) - Av”^ a+1) T _(a+1) T+B(v+T)+Cv+A-2B-C = -^Xx 2 , 
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(7.11) - Av _(a+1) x~ (a+1) (v+x)+(2 B+C) (v+t)+2A- 

-2 (2B+C) = -Avx . 

We now solve (7.11) for -A(vt) and substitute the 

resulting expression into (7.10), obtaining an equation with 
a factor v - t . Since we wish to study the locations of 
possible solutions of (7.1), (7.2) with v ^ t , we 

cancel this factor obtaining 

(7.12) (B+C) (v+T) + A -2B - C = % Ax 2 . 

Let us denote the solution of (7.12) for v in terms of 
t by v#(x,A) . If we replace v in (7.10) with v#(x) , 
we obtain 

(7.13) Ax" (0t+1) = (Ax-C) [ v # (x,A)] a+1 

For x > 0 , trivial solutions of our boundary value 
problem are characterized by v = x = M A) £ (0,1) , 

N = T =- 1 sxk(x) 2 .We accordingly limit our attention to a 
neighborhood (possibly quite large) of the trivial solutions 
by seeking solutions of (7.13) for which 

(7. 14a,b,c,d) x e (0,1), v* (x,A) e (0,1) , 

N(v# (x , A) , x) < 0 , T(v*(x,A),x) < 0 

for A > 0 . We wish to characterize where (7.13) has 

solutions satisfying (7.14) . Such solutions lie in the 
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complement of G . (If (7.14) is not imposed, then (7.13) 
always has a solution.) 

Note that (7.14c) is automatically satisfied by 
virtue of (7.1) . We now use (7.12) — (7.14) to show that 
the set of points (v,t) e (0,1) x (0,1) at which (7.1) , 
(7.2) has solutions for (7.8), (7.9) can be quite small. 

Call this set H(X) . Clearly 

(7.15) H (X) = { ( v , t) £ (0,1) x (0,1) : v = v # (x,A) , 

x satisfies (7.13), T(v#(x,A),x) < 0 } 

We now obt'ain a number of inequalities that elements of 
H( A) must satisfy . 

Equation (7.3) implies that (7.14d) is satisfied 
if and only if 

(7.16) v # (t , A) > t/ 2 
Equation (7.12) implies that 

(7.17) v # (x , A) > -x+ 1 + § -•-= -§ 

which can represent a very severe or an insignificant 
restriction depending on the relative size of A, B, C. 
Equation (7.12) and conditions (7.14a,b) imply that 

x < [ 2 (C+A) ] ^ A -J * 


(7.18) 
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A simple graphical analysis of (7.13) based on (7.14a,b) 
shows that 


(7.19) 


r 

x > min 


A+C 

X 



Note that if B-A and X are sufficiently large, 
then the set of t's satisfying (7.14), (7.17), (7.18) is 
empty. If X is small enough, then (7.19) implies that 

X >_ 1 

By substituting (7.16) into (7.12) we obtain 


(7.20) x t [T_,T + ] , Xt + = |(B+C)±[f (B+C) Z +2X(A-2B-C)] J5 


the interval [x_, x ] being regarded as empty if x + are 

not real. If A-2B-C > 0 and if X is large enough, then 
(7.20) delivers a lower bound for x that is larger than 
that of (7.19) . 

We seek similar estimates independent of X . Let 
us use (7.18) to get a lower bound for Ax . a 

sketch shows that a lower bound for x is the solution of 


(7.21) 


A [ 


2 (C+A) - (ot+1) / 2 
X J 


(Xx-C) 


Thus 
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(7.22) 


T > 


AX 


(ot-D/2 


= X 


[ 2 (C+A) ] 


(a+l)/2 


If a e (0,1), then (7.22) delivers a lower bound depen- 
dent on X , but nevertheless better than those of (7.19) 
and (7.20) when X is large enough .If a = 1 , then 

(7.22) gives bound that is clearly better than (7.19) and 
that itself has a further lower bound independent of X . If 
a > 1 , then 

2 t > £x -2 + (l-£)x a_1 

(7.23) 

= ( a+ l )C (ot - 1)/(a+1) A 2/(a+1) 

(a-1 ) (o-D/(a+l) 2 -2(a+2)/(a+l) {C+A) 

with x = [ ( a 3i) A 1 ^ ^ a+1 ^ • C = C + A > the expressions 

on the right side of (7.23) being the minimum of the right 
side of (7.22) for X e (0,») . An upper bound for the 
rightmost term of (7.23) can be obtained by the Young 
inequality. The right side of (7.23) can be made arbitrarily 
large by taking a -1 sufficiently small and positive. 

It is not difficult to obtain further estimates, e.g., 
estimates for the neighborhood discussed in the comments 
following (7.7b). We refrain from doing so. 

The basic message of this section is that the set G 
can be very large and that precise information about it can 
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be gained directly from degenerate constitutive equations. 
Moreover, we would be able to strengthen the global results 
of Section 7 even more by obtaining estimates for the 
location of solutions of our boundary value problem and res- 
tricting G to refect such estimates. We do not attempt such 
estimates here. 
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8. CONCLUSION . 

In an unshearable shell n is constrained to be zero. 
The shear resultant H , not defined constitutively , is the 
Lagrange multiplier maintaining the constraint of unshear- 
ability. It is a fundamental unknown of the constrained 
problem. The same analysis as that used above can treat 
unshearable shells if H is treated as n 1S * 

Some of our techniques are reminiscent of those of 
Antman (1978) for plates. The difference in boundary con- 
ditions forced us to give the coordinate functions r and z 
a central role in the analysis of Section 5, which they 
lacked m the plate theory. The issues treated m Section 
7 never arose for plates. Of course the nodal properties 
of branches for shells differ markedly from those for plates 
and for circular arches (cf. Antman & Dunn, (1980).) Incident- 
ally, the uniqueness theorem of Antman (1978, Section 7) is 
faulty; it can be readily corrected by following the develop- 
ment of Section 6 above. 

It is illuminating to compare our results with the 
numerical results of Bauer, Reiss, & Keller (1970). Treating 
an engineering theory of shells they discovered the remark- 
able fact that all the bifurcating branches they analyzed 
are connected. Our results point in the diametrically 
opposite direction. We showed that nodal properties of 
branches bifurcating from simple eigenvalues could only lose 
their nodal structure at certain special places in solution- 
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parameter space. In particular, if we could obtain delicate 
pointwise estimates showing that v c (0,1), t e (0,1) , 

N < 0 , T < 0 , everywhere on a family of bifurcating 
branches, then the results of Sections 6 and 7 show 
that nodal properties cannot switch along such branches. 
Consequently such branches could not connect bifurcation 
points with simple eigenvalues having different eigenfunctions. 

Now Section 7 shows that for fixed X the comple- 
ment of G is generically a finite number of points. On 
the basis of such genericity arguments we could expect that 
on a given bifurcating branch there would be no s 0 e [0,ir] 

and X >0 such that (v(so) , t(So)#X) t G . ( Such generic- 
ity arguments are dangerous: The nodal structure we found 
for n and \p-s is not generic for pairs of functions. ) 

The chief source of the suspected variance between 
our results and those of Bauer, Reiss, & Keller is likely 
the fact that various geometric terms regarded as small 
were discarded in their theory ( as is typical of all tech- 
nical shell theories) . Consequently these terms could not 
come into play in the large solutions they computed. A 

conceivable alternative explanation of the connectedness 
of all the bifurcating branches found by Bauer, Reiss, & 

Keller is that they correspond to spurious solution branches 
of the sort often found in numerical studies of bifurcation 
problems. We tend to dismiss this explanation because it 
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does not seem compatible with the consistent pattern they 
found. Bauer, Reiss, & Keller did not qive enouah data for 
us to determine what sort of nodal properties their solutions 
possess. 

Our approach may also be compared with that of 
Lange & Kriegsmann (1981) who studied continuation prop- 
erties for very thin shells by using asymptotic expansions. 
Our global bifurcation theory gives an approach to the con- 
tinuation of branches that does not require special hypoth- 
eses on thickness, or more generally, on constitutive func- 
tions. Presumably their detailed approach could be applied 
to our geometrically exact theory for constitutive functions 
corresponding to very thin shells. We could also introduce 
the thickness h as a second bifurcation parameter and 

study the role of thickness by using global multiparameter 
bifurcation theory (cf. Alexander & Antman (1981)). 

We mention that the perturbation results of 
Shih (1985) show that the local disposition of bifurcating 
branches depends critically on the nature of our nonlinear 
constitutive equations. In particular, the local structure 
can be quite different from that found in technical shell 
theories, which essentially use linear constitutive relations 

As we have indicated above, the process of numerically 
constructing solution branches can be quite tricky. It would 
seem that nodal properties on solution branches could give 
an extra control on numerical computation. But, to our 
knowledge, no such use of nodal properties has been attempted 
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§ 9. APPENDIX A. FORMULATION OF THE GOVERNING EQUATIONS . 

In this appendix we describe a special director theory 
for a nonlinearly elastic axisymmetric shell undergoing an 
axisymmetric deformation. Let {i,j,k} be a right-handed 
orthonormal basis for the Euclidean space E 3 

In the special Cosserat theory we employ, a shell is 
regarded as a surface together with a unit vector field, 
called the director field, defined on it. For such an 
axisymmetric shell undergoing an axisymmetric deformation, the 
surface is generated by a plane curve rotated about the 
k-axis and the director field is likewise generated by 
rotating a director field coplanar with the generating curve. 



Fig. 9-1 
Reference Configuration 


Let 


(9.1) 

(9.2) 


ei($) = cos <j>i + sin <pj 
e 2 ( <J) ) = -sin <J)i + cos <pj 
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The position r*(s,<j)) of a typical material point 
on the (mid) surface of an axisymmetric shell in its unde- 
formed reference configuration is (see Fig. 9-1), 

(9.3) r*(s,4>) = r*(s)ei(<j>) + z*(s)k , 

where <f> is the meridional angle and s measures arc length 
parameter. (s,<p) identify each material point of the 
(mid) surf ace . We let s e [0,L] . Define 

(9.4) a*(s,<f>) = cos ip*(s)ei(<p ) + sin ^*(s)k , 

(9.5) b*(s,cf>) = - sin \p*(s)ei(<p) + cos ^*(s)k 

where ^*(s) is the angle between the direction of the tan- 
gent vector to the curve r*(*,0) at s and that of ei(<j>) . 

Deformed Configuration 

The axisymmetrically deformed configuration of an 
axisymmetric shell is determined by a pair of continuously 
differentiable functions 

(9.6) [0 ,L] x [ 0 , 2 tt] s (s,4>) r (s, <J>) , b(s,<J>) 
where 

(9.7) r(s,4>) = r(s)ei(4>) + z(s)k , 

(9.8) b(s,<Ji) = -sin (s) ej (<J>) + cos i^(s)k 
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Define 

(9.9) a(s,tj>) = cos ij; (s) ej. (<J>) + sin ^(s)k , 

(9.10) r , (s, <j>) = v (s) a (s, <£) + n(s)b(s,<|>) , 

~ O V -V 

(9.11) r , ^ (s, <f>) = r(s)e 2 (<j>) , 

where r, (s,<J>) and r, ,(s,<j>) are partial derivatives of 

~ b ~ (p 

r(s,<J>) with respect to s and <p . (See Fig. 9-2) 


Fig. 9-2. 
Deformed Configuration 

The strain variables ( which do not all vanish in the 
reference configuration ) consist of invariants ( namely the 
scalar product and scalar triple products ) of 

r , (s, <}>) = va (s, tj>) + nb(s,<j>) ; 

O *V -V 

b, (s,<j>) = — xjj ' a ( s , <J) ) ; 



(9.12) 
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r , ^ (s, <(>) /r* (s) = [r (s) /r* (s) ] e 2 (<£) ; 

b, ^ (s, <j>) /r* (s) = - [siniHs) /r* (s) ] e 2 {<(>) 

These invariants depend only on 

• (9.13) w ( s ) = (v,x,n,y,a) , 

where 

(9.14) T (s) = r (s) /r* (s) ; 

(9.15) y (s) = ip' - (t^*) ' ; 

(9.16) a(s) = sin ifi(s)/r*(s) - sin if>* (s) /r* (s) 

and the prime denotes the derivative with respect to s 

REMARKS . (1) In the case of a spherical shell of radius 1, 

where 4**(s) = s , r*(s) = sin s , we have v = t = 1 , 

11 = 0,0 = y = 0 m the reference configuration. Thus 

v - 1 is a measure of longitudinal strain, x - 1 the 
meridional strain, n the shear strain and a and y are 
the circumferential and longitudinal bending strains respec- 
tively. Note that y + 1 and o + 1 are almost curvatures. 

(2) In the uniform contraction and expansion of 
a spherical shell into one of different radius, v - 1 and 
x - 1 change but y and o remain zero, so that y and 
a isolate the effects of bending. 
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EQUILIBRIUM EQUATIONS . 

We suppose the shell is subjected to a uniform external 
hydrostatic pressure of constant intensity x per unit actual 
area. Accordingly, the corresponding force per unit reference 
area of the shell is 


(9.17) 


X 






r x 
~,s 





X 


r x 



r* 


Let 


ni be contact force resultant per unit reference length 
acting across s = constant. 

n 2 be contact force resultant per unit reference length 
acting across <p = constant. 

be contact couple resultant per unit reference length 
acting across s = constant. 

m 2 be contact couple resultant per unit reference length 
acting across tp = constant. 

The equilibrium of total force on the material of 
{(s,<j>): s x < s < s 2 ; $i < <t> < <l >2 } implies 


(9.18) 


'<t>2 


<t> 1 


[n 1 (s 2 ,<j>)r*(s 2 )-n 1 (s 1 ,(J>)r*(s 1 )]d<t> + 


f sz 

[n 2 

Si 


'02 

'S 2 


Si 


5-s x £•* 

X r* d4>ds = 0 
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If n 2 and n 2 are continuously differentiable, we can 
differentiate (9.18) with respect to <j> 2 and s 2 and then 
replace <J> 2 and s 2 with <p and s to get 

Sn 2 

(9.19) [r* (s)ni ( s , 4> ) ] , g + (s,<j>) +Xr, g (s,<j>) * 


x r , ( s , 4 > ) = 0 


The equilibrium of total moments on the same shell element 
implies 



(<t>2 

s 2 

(9.20) 

(r x m + mi ) r* (s) 

d<J> + 

• 

*1 ~ ~ 

Sl 


rs 2 


(r x n 2 + m 2 ) 


ds + X 


<J>i 


<t>2 

<Pl J 


r s 2 


Sl 


Si 


r x (r, s x r,^)dsd<}> = 0 


Thus, if mi and m 2 are continuously differentiable, we 
obtain 


(9.21) 


[r x mr*(s)], s + [r*(s)mi], s + (£ x n 2 + m 2 ),^ + 
+ Xr x (r, g x r,^) = 0 


Substituting (9.19) into (9.21) , we get 


(9.22) 


[r* (sjmo. (s,4>) ] , s + r, g (s,<J>) x r* (s) nj. (s, <j>) + 
+ m 2/(( )= 0 
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We seek axi symmetric solutions, for which we require 


(9.23) 

ni (s , 4> ) 

/A 

= N (s) a (s, 4>) + 

A 

H (s) b (s , <j>) 

/ 

(9.24) 

n 2 ( s , <f> ) 

A 

= T (s) e 2 (<j>) 

r 


(9.25) 

mx ( s , 4> ) 

/A 

= - M ( s) e 2 (<j>) 

r 


(9.26) 

m 2 (s, <f>) 

• e 2 = 0 

• 



Thus 

m 2 = (m 2 • a) a + (m 2 • b)b with components regarded 
independent of <j> . Then 

9m2 

(9.27) —— = (m 2 • a)cos ^e 2 ( ) - (m 2 • b) sin \pe z 

a (p ~ ^ -v ~ ~ - 

= k x m 2 

Thus (9.19) and (9.22) become 

(9.19a) {r* (s) [N (s) a (s, <j>) +H (s)b(s,<f>) ] } , s - T(s)e 2 (4>) + 

+ Ar(s)r(s,<j>) , x e 2 (<j)) = 0 

~ O ~ 'V 

(9.22a) [r* (s)M(s) ] f _+r* (s) [vH(s) -nN(s) ] - (k x m 2 ) • e 

S -v «w 

A 

Now we set (m 2 • a)cos ^(s)-(m 2 * b) sin \p(s)= Z (s)cos ip(s) 
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to obtain 


(9.22b) 


A A A A 

[r* (s)M(s) ] '+r* (s) (vH-nN) - E(s)cos iMs) = 0 


SPECIALIZATION TO SPHERICAL SHELLS . 

Without loss of generality we take the shell radius 
to be 1. Thus r*(s) = sin s, ip * ( s ) = s , z*(s) = -cos s, 

L = 7T , and equilibrium equations (9.19a) and (9.22b) 
become 

A A /V 

(9.28a) [N(s)sin s]'-T(s)cos iMs)- Xrri -H(s)sin s^'(s)=0 , 

A A 

(9.28b) [H(s)sin s]'+T(s)sin i|j(s)+ Xrv +N(s)sin s^'(s) = 0 , 

A A A /V 

(9.28c) [M(s)sin s}'-E(s)cos ip (s) + ( vH-r|N) sin s = 0 

In order to understand the component Z occuring in 
(9.22b) , it is convenient to have recourse to the principle 
of virtual work. 


Let y(s,<f>), x(s,<J>) e C 2 (ft) with compact support, 

where Q = [si,b] x [<}>!, 6], 0 < b < L, 0 S ^ 2 tt . The 

principle of virtual work requires that 


(9.29) 0 = 




Sl 


{ r*n x *y , s +n 2 *y , (rxnj+mz ) r*x, g + 


+ (rxn 2 +m 2 ) x, -X (r . xr , ) *y-Xrx (r , _xr , ) .x}dsd<{> 
~ ~ - m ~ a ~ s ~ <P ~ ~ ~ s 
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It is convenient to use the traditional notations. Let 
Sv , 6 n and 6\p be in C 2 (ft) with compact supports and 
set 

(9.30) y = Sr - (r*e 2 )SiJ) ; x = -S^e 2 , 

where 

(9.31) (Sr) , = Sr. = Sva + Snb + (vb-na) Si|) 

•v b ~ b *S» *v ~ ~ 

Thus , 

(9.32) Sr,^ <5re2 = (Svcos ii>-6nsin i|;-vsin iJ;5^-ncos ipSip)e 2 
Also, we get 

. _ Sr 

(9.33) 0T r* 

(9.34) 5a = ~^~jr cos ip , 

r 

(9.35) 6\i = S\p • , 

(9.36) X '<p = 


Thus (9.29) becomes 
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( 9 . 37 ) 


0 


' 6 rb 

(Ni5v+H<5n+T6T+M6y+I6a) r*dsd<£ - 

J Si 


-X 


r6 b 
4>i-' s ; 


(r , g - x r,^) • 6r dsd$ 


This is the traditional form of the principle of virtual work 
with E as the generalized force corresponding to the strain 
variable a 


Note that in the expression (m 2 • a) cos ip- (m 2 • b) 

A /\ 

sin \p= Ecos ip, we do not imply that E = m 2 • a and that 


b • m 2 


0 
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